Draft version February 1, 2008 

Preprint typeset using lAT^rjK style emulateapj v. 04/03/99 



o 
o 
o 

(N 



COOL COMPANIONS TO HOT WHITE DWARFS 
Paul J. GreenQ 

"^Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Science Foundation. 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 

email: pgreen@cfa.harvard.edu 



Babar Ali^ 

Infrared Processing and Analysis Center, Caltech, Mail Code 100-22, 770 S. Wilson Avenue, 

Pasadena, CA 91125 
email: babar@ipac.caltech.edu 

p |. AND 

I R. Napiwotzki 

■ Dr. Remeis-Sternwarte, Sternwartstr. 7, 96049 Bamberg, Germany 

I email: ai23@sternwarte.uni-erlangen.de 

^ I Draft version February 1, 2008 

O : ABSTRACT 

I . Low mass companions to high mass stars are difficult to detect, which is partly why 

'nI" \ the binary fraction for high mass stars is still poorly constrained. Low mass companions 

^ ' can be detected more easily however, once high mass stars turn into white dwarfs. These 

O . systems are also interesting as the progenitors of a variety of intensely-studied interacting 

binary systems, like novae, CVs, symbiotics, Ba and CH giants, Feige 24-type systems, 
Ph' and dwarf carbon stars. 

Q . We describe a near-IR photometric search for cool red dwarf companions to hot white 

^ \ dwarfs (WDs). IR photometry offers a sensitive test for low mass main sequence (MS) 

^ ■ companions. Our sample of EUV-detected WDs offers several advantages over pre- 

vious (largely proper motion-selected) WD samples: (1) the high WD temperatures 
; (24000 < Teff < 70000K) insure excellent IR flux contrast with cool dwarfs (2) the 

^ ■ range of evolutionary parameter space occupied by the WDs is considerably narrowed 

and (3) the random effects of the intervening ISM provide a complete but reasonably-sized 
sample. While some composite systems have been found optically among WDs detected 
in recent EUV All-Sky Surveys, we develop an IR technique that probes farther down 
the main sequence, detecting yet more companions. 

We use detailed DA model atmosphere fits to optical spectra to predict K magnitudes 
and distances, against which we contrast our near-IR observations. Our photometric sur- 
vey reveals 10 DAs with a significant excess in both J and K. Half are newly discovered, 
and are most likely previously unrecognized binary systems. Neither the frequency of 
infrared excess nor the mass estimate of the red dwarf companion correlate with white 
dwarf mass, as might be expected if either the EUV detectability or mass of the white 
dwarfs were significantly affected by a companion. Infrared spectra of these systems 
should help to determine the mass and spectral type of the cool companions presum- 
ably causing the IR excess, leading to better estimates of the mass ratio distribution in 
binaries. 

Counting previously known binaries, and resolved pairs, we find the total binary fraction 
of the sample is at least a third. Since most WD progenitors had initial masses > 2Mq, 
we thus provide a photometric measure of the binary fraction of high mass stars that 
would be difficult to perform in high mass main sequence stars. We estimate that 90% 
of the companions are of type K or later. 
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1. INTRODUCTION 

EUV-detected stars have revealed in their UV 
spectra the presence of about 15 hot white 
dwarf (WD) companions to bright stars in 
non-interacting binary systems (e.g., Burleigh, 
Barstow, & Fleming 1997). These WDs are hid- 
den at optical wavelengths, because of their close 
proximity to much more luminous companions - 
main sequence stars of spectral type K or earlier, 
or evolved stars. 

Many interacting binary systems where the WD 
is the primary (i.e., optically brightest) star have 
also been found among EUV-detected systems 
(e.g., 6 close, interacting white dwarf/red dwarf 
binaries by Vennes & Thorstensen 1994). Opti- 
cal or ultraviolet spectral observations are most 
commonly used to detect companions to WD pri- 
maries, by searching for (1) the presence of narrow 
Balmer line emission overlying the broad smooth 
Balmer absorption of the WD, (2) a composite 
WD -|- main sequence spectrum, or (3) radial ve- 
locity (RV) variations. However, only WDs with 
very close, or intrinsically active companions will 
be found by method (1). For hot WD systems, 
composite spectra (2) are only expected to be 
visible if the companion's spectral type is early. 
RV variations (3) require multiple observations at 
high spectral resolution, and detection strongly 
favors close and/or massive companions. 

Such discoveries have been strongly dominated 
by these selection effects, with companions bi- 
ased to earlier types than predicted by the sim- 
ulations of deKool & Ritter (1993) and others. 
Scaling from deKool & Ritter's results, Vennes & 
Thorstensen (1994) estimate that "at least twice 
as many close binary systems remain to be identi- 
fied from EUV surveys, most of them with a low 
mass secondary." The resulting sample of bina- 
ries known to date therefore must diverge strongly 
from the intrinsic distribution, in overall normal- 
ization, as well as in mass and spectral type of the 
main sequence companions. 

Large proper motion-selected samples (includ- 
ing those observed in the IR e.g., Probst 1983a,b 
and Zuckerman & Bccklin 1992), are kinemati- 
cally biased against luminous (statistically more 
distant) young, hot WDs, and span a bewilder- 
ing range of temperatures, ages, and metallici- 
ties. Our EUV-selected sample is instead dom- 



inated by freshly-minted disk DAs, providing a 

narrower slice of parameter space for stellar evo- 
lution (in the wide systems) or for dynamical evo- 
lution (in any new close systems). Objects known 
or proposed to contain white dwarf (WD) stars 
in currently or previously interacting binary sys- 
tems constitute diverse laboratories for the study 
of stellar evolution. A partial hst includes novae, 
cataclysmic variables, symbiotic stars, Ba and CH 
giants, Feige 24-type systems and dwarf carbon 
stars (Green & Margon 1994). 

We describe here a near-IR photometric sur- 
vey for low mass companions to hot white dwarfs 
(WDs). By investigating only EUV-detected 
WDs, we obtain a very reasonably-sized but com- 
plete sample of young WDs, next to which very 
late-type dwarf companions can be detected in the 
near-infrared by searching for a K excess. Many 
hot white dwarfs {T^s > 24000K; Finley et al. 
1993) have been detected in the recent EUV all- 
sky surveys. EUV detection of these hot WDs 
depends primarily on their temperature, distance, 
and the intervening Galactic ISM. Our sample of 
EUV WDs (whose selection we define below) of- 
fers excellent flux contrast in the IR relative to 
optical; cool companions will almost always be 
brighter in the K band than the hot WDs. IR 
photometry of hot white dwarfs is nearly non- 
existent in the literature. 

To know what K mag to expect for the WDs, we 
benefit from constraints on logg, radius, and Teg 
derivable from optical spectra for the WDs in our 
sample, using NLTE model atmosphere fits (Napi- 
wotzki. Green, & Saffer 1999; hereafter NGS99). 
The resulting predictions for K magnitudes allow 
a direct search for any IR excess from a cool com- 
panion. In some cases, IR colors also provide a 
preliminary spectral type for the companion. 



2. SAMPLE SELECTION CRITERIA 

We chose to limit our initial sample to DAs, for 
which model atmospheres provide the best tem- 
perature and mass constraints. We start with 
73 known DA white dwarfs in the EUVE bright 
source list (Malina et al. 1994). We exclude 
sources at low galactic latitudes (|6| < 15) to 
avoid crowding problems, and in the South {5 < 
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—20), yielding a list of 27 DAs. ^ A similar pro- 
cedure for DAs listed in the i?0S'.4rWide Field 
Camera survey Bright Source Catalogue (Pounds 
et al. 1993) culls an additional 30 objects. The 
Table [| lists the full sample of 57, including pub- 
lished data used in this paper. Model atmosphere 
fit parameters for DAs that are not available in 
NGS99 are taken from Vennes et al. (1997a) or 
other references listed in Table |l]. 

Of the 57 sample objects, we present new 
IR photometry for 47. We did not ob- 
tain new measurements for 10 stars with pub- 
lished sensitive optical spectrophotometry and 
IR photometry and/or known binaries from 
optical studies (Feige 24, HZ 43, CD 50, 
V471 Tau (Vennes, Christian, & Thorstensen 
1998), PC 0824+289 (Heber et al. 1993), 
HD74389B (Liebert, Bergeron, & Saffer 1990), 
REJ 1016-052, REJ 1426+500 (V97), and IK Peg 
(Wonnacott et al. 1993). Due to observing con- 
straints (a combination of weather, poor seeing 
and faint objects, or celestial placement of ob- 
jects) no photometry was obtained for 0138+252 
(PC 0136+251). 

On the other hand, we report new IR data 
for 4 objects that fell outside the strict sam- 
ple definitions just outlined. These include the 
known binaries REJ 1036+460 (PC 1033+464), 
REJ 1629+780 as well as REJ 0148-253 
(CD 1401, also a known binary), and REJ 0457- 
280 (MCT 0455-2812) which are outside the sam- 
ple declination limits. 

We note that since the initial sample selec- 
tion, several relevant discoveries pertaining to 
sample objects have been made. REJ 0134-160 
(GD984) has central Balmer emission components 
produced by a dMe companion. RE 1440+750 
turns out to be a magnetic DA (Dreizler et al. 
1994). 

3. OBSERVATIONS AND DATA REDUCTION 

We obtained J- and K-haiad photometry of the 
sample objects at the Kitt Peak National Obser- 
vatory^ 2.1-m telescope using the InfraRed IM- 
ager (IRIM; Fowler et al. 1988). The 256x256 
IRIM, completed in 1992, is a HgCdTe NIC- 
M0S3 array, which at the //15 focus of the 2.1- 
m yields a 280" x 280" field of view (plate scale of 

-'^ Emission from 40 Eri was resolved with the Einstein HRI, and likely arise mostly from 40 Eri C, the dMe flare star 
(Cash et al., 1980), not the white dwarf 40 Eri B. 

^The Kitt Peak National Observatory is operated by the Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National Science Foundation. 



1.09"pixel-i). 

Table |^ & ^ summarize our observing proce- 
dure. Table lists the individual observing runs 
dedicated to this project. Column (1) identifies 
the run. The Universal Time (UT) date of the 
run and the photometric standard stars used dur- 
ing each night are shown in Columns (2) and (3). 
Column (4) identifies the references containing in- 
frared (J- and K- band) magnitudes for the stan- 
dard stars, and column (5) contains a short over- 
all subjective description for the night conditions. 
Table ^ summarizes our observing strategy. Col- 
umn (1) identifies the program object (given as 
an RE J nnnn+mmm) . Column (2) lists the nights 
on which the target was observed, and column (3) 
any exceptions to our standard observing proce- 
dure of 6 dithered images per target, each image 
consisting of 6 co-added 10 second exposures. The 
exceptions are noted in column (3) of Table ^. 
The number of exposures depended on the faint- 
ness of the object. To minimize any hysteresis ef- 
fects in pixel sensitivity, during the observations, 
we varied the direction and distance between the 
position offsets, and never placed a target at the 
center of the array. At least once during each 
night, 10 - 20 (co-added) dark frames of 10 sec. 
duration were also obtained. These are median- 
combined to form a single dark frame free of cos- 
mic ray hits. The atmospheric seeing (and hence 
the FWHM of stars) varied from 1.4-2.5 pixels 
during the observations. 

We reduced our image data using the Image Re- 
duction and Analysis Facility (IRAF) software. 
For each filter separately, the reduction proceeds 
as follows. All images are corrected for any de- 
tector non-linearities (R. Joyce, private commu- 
nication). Next, any bad pixels (as defined by the 
bad pixel mask; see below) are replaced with the 
average of the linearly interpolated values from 
adjacent pixels in both dimensions of the image. 
Each image is then sky subtracted and divided by 
the fiat-field. 

The sky frame is produced by median- 
combining all exposures of the star. This also 
serves to subtract the dark current, since the dark 
current is included in the sky frames. The fiat- 
field image is created as a median-combined frame 
from a large (> 30) number of individual images. 
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Individual images are median-subtracted to ac- 
count for sky/background variations. The flat- 
field image is dark subtracted and normalized. 
The flatfield image is also used to produce the 
bad pixel mask, as defined by the histogram dis- 
tribution of pixel intensities. All pixels with in- 
tensity values near or beyond the extreme values 
of the histogram distribution are defined to be 
bad. The bad pixel mask was produced once and 
used throughout. 

Since our sample is defined to be at high Galac- 
tic latitude and hence uncrowded, simple aper- 
ture photometry is adequate. Following the ar- 
guments presented in Howell (1989) we used 
a 4.0 pixel diameter aperture to maximize the 
signal-to- noise ratio (S/N). Instrumental magni- 
tudes are measured for each object using the 
IRAF task QPHOT with the radius of the inner 
sky annulus set to 10 pixels and the width of sky 
annulus set to 5 pixels. 

For target stars REJ 0443-034, REJ 0907+505, 
REJ 0940+502, and REJ 1442+632, visual inspec- 
tion failed to reveal the object on individual im- 
ages. In such cases, we shift (by integer pix- 
els) and combine all the images of the field to 
form an average mosaiced image. Individual 
frames are median-subtracted prior to averaging 
to remove any zero point offsets due to varying 
sky/background conditions. The pixel shifts are 
calculated from other nearby stars in the field. 
The photometry procedure for these mosaiced im- 
ages is described in Appendix ^ 

We derive the zero point calibrations for the 
instrumental magnitudes from photometric stan- 
dard stars observed at a variety of airmasses each 
night. These standards stars were mostly selected 
from the UKIRT faint standard list (Casali & 
Hawarden 1992). The photometric standards for 
each night are listed in Table ^ 

We eliminated the most extreme outlying mea- 
surement for an object if its elimination reveals it 
to be > 5x the root-mean-square (RMS) of the 
resulting average magnitude. The resulting rela- 
tive photometric error a is taken to be the RMS of 
all measurements from each individual frame. To 
obtain the final magnitude error, we add to a the 
RMS dispersion in the standard star zeropoints 
taken over the night. The final magnitudes for 



targets are the mean values produced by averag- 
ing the photometry from all nights on which the 
target was observed. If the target was observed 
on more than two nights, we use the algorithm 
described above to eliminate outlying measure- 
ments. 

Table ^ lists for each object with new IR ob- 
servations both the observed and predicted mag- 
nitudes. Column (1) lists the RE name, (2) the 
observed J magnitude and (3) error. Columns 
(4) and (5) list the K magnitude and error, re- 
spectively, with the (J — K) color and error in 
columns (6) and (7). The absolute magnitude 
My from the best-fit spectral model is shown in 
column (8). We list in column (9) the predicted 
{y — K) color based on synthetic photometry of 
the best-fit DA models. Predicted near-IR mag- 
nitudes for the DA WDs are listed in columns 
(lO)-(ll), and notes in column (12). We note that 
our predictions of near-IR colors for these systems 
have very little scatter over the full sample range, 
with mean (J — K) = —0.24 ± 0.02, as can be 
seen from Figure |I[ DA colors - even {V — K) - 
are nearly degenerate for DA temperatures higher 
than about 35,000K. 

4. ANALYSIS 

4.1. Predicting Near Infrared Magnitudes and 

Colors 

We use a grid of 184 NLTE model spectra for 
DA WDs, with T^g from 20,000 to 100,000 (in 
16 steps, coarser toward the the high end) and 
logg in 12 equal steps from 7.0 to 9.75. We pre- 
dict colors or magnitudes for each grid model by 
convolving the model spectrum over the observed 
filter bandpasses. In the visual, we use Mathews 
& Sandage V filter curves and typical CCD quan- 
tum efficiency (QE) curves, while in the near-IR, 
we use the actual NOAO filter transmission curves 
for the J and K filters, along with the IRIM 
QE curve. Zeropoints for all magnitudes are cali- 
brated by convolution with the Vega spectrum of 
Hayes & Latham (1975), renormalized (by a fac- 
tor of 1.028) to zero magnitude. As a check on our 
synthetic photometry, we tabulated synthetic col- 
ors derived from the spectrophotometric atlas of 
Bruzual, Persson, Gunn, & Stryker (BPGS; Gunn 
& Stryker 1983) which includes infrared (IR) dis- 
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Fig. 1. — Sample DA white dwarf temperature from best- fit model atmospheres of NGS99 plotted against our pre- 
dicted (J — K) color for the star assuming it is solitary. DA colors are nearly degenerate for temperatures above about 
35,000K. 



tributions from Strecker et al (1979). Our syn- 
thetic colors agree with published colors for the 
corresponding spectral types to better than about 
0.05mag for blue spectral types (Zombeck 1990). 

We predict colors for each program WD star 
from its best-fit gravity and temperature using 
a 2-dimensional linear interpolation between the 
colors derived for the nearest grid models. Pre- 
dicted colors for each program star are combined 
with observed visual DA magnitudes, to predict 
their J and K magnitudes, assuming that they 
are isolated. A table of Jp and Kp (predicted J 
and K magnitudes) for DA stars are listed in Ta- 
ble ^ along with the observed magnitudes. 



4.2. Deriving Distances and Searching for IR 

Excess 

Because each of the model spectra are pre- 
sented in stellar surface fluxes Fx (in cgs units 
erg cm~^ s~^ cm~^), magnitudes we derive for 
each model in the NLTE grid are effectively un- 
normahzed surface magnitudes Vs = 2.5 log F^. 
The formula we derive for converting the best fit 
DA WD surface fluxes Fx (in erg cm~^ s~^cm~^), 
surface gravities log (7, and masses into abso- 
lute magnitude is 



Mv = -2.51ogFA+2.51ogc/-2.51og +31.042 

Our own resulting values of My differ little from 
those for 35 DAs in common with Vennes et al. 
1997a (mean difference 0.01 ± 0.24mag). Since 
My correlates strongly with logg, this is consis- 
tent with the more detailed comparison performed 
by NGS99 of internal, external, and systematic er- 
rors in logg and Teg between the two studies. For 
distances to the DA WD, we compare to tabulated 
V magnitudes, mostly from Table 3 of Marsh et 
al. (1997). 

We derive the i^' magnitude of any cool compan- 
ion, Kc, from the difference between the predicted 
magnitude of the DA WD, Kp, and the total ob- 
served magnitude Ko using 

Kc = Kp- 2.5 log [dexO A{Kp - Ko) - 1] 

In Table |] we list the properties of candidate 
cool companions detected in our IR-observed sam- 
ple of DA WDs, estimated using the combination 
of the best-fit NLTE models of NGS99, our own 
infrared photometry, and published optical mag- 
nitudes. Column (1) lists REJ name. The com- 
panion J and K magnitudes calculated as detailed 
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above are listed in columns (2) and (3). Spectral 
type of the candidate companions can be crudely 
estimated from the derived M^^ or from J — K. 
Since the combined errors on M^^ are in general 
lower, we list this in column (4) and the resulting 
spectral type estimates in columns (5) and (6), as- 
suming young or old disk population for the red 
dwarf. The implied V magnitude (from V — K 
for the old disk spectral type) is shown in column 
(7), with notes in column (8). 

4.3. Candidate Cool Companions 

Our adopted criteria for claiming a significant 
near-IR excess over the predicted values for an iso- 
lated DA is simply a 3cr excess in both J and K 
relative to the predicted IR magnitudes of the DA 
white dwarf, Jp and Kp. Errors in our observed 
K magnitudes are generally 2-3 times higher than 
for J, but provide greater wavelength contrast to 
the optical. We do not use {J — K) as a dis- 
criminant because its combined errors are highest. 
The errors we adopt for the excesses {Jp — Jo) and 
(Kp — Ko) include our photometric errors in J and 
K, combined with assumed errors of O.lmag each 
on the predicted values of Kp and Jp. 

As an example, for REJ 1845+682, we calcu- 
late from the published V magnitude and model 
atmosphere fits in Table 1 a predicted magnitude 
with its estimated error of Kp = 16.2 ± 0.1 (Ta- 
ble 4). The difference between this and the ob- 
served magnitude K = 14.37 ± 0.19 (Table 4) is 

thus Kgxcess 

8cr excess. 



1.83 ± 0.21mag, qualifying as an 



Of the 47 objects with reliable K photometry 
(which excludes REJ 1235+233/PG 1232+238) we 
find 10 systems with a significant IR excess. Of 
these, 5 are known binary systems. Of the 5 that 
are newly recognized systems with IR excess, 4 
have spatially resolved objects discernible either 
in the near-IR or optical images within 10". In 
only some cases, noted in Table ||, these objects 
are included in the near-IR photometry. 

The overall properties of DAs with IR excess 
reveal no obvious bias within our observed sam- 
ple. The mean distance of the DAs without ex- 
cess is 144±98pc while those 11 with excess have 
117 ± 53pc. Similarly, the mean values and dis- 
tributions of Teff , log (7, V, and My of the DAs 
without excess are indistinguishable from those 
11 with excess. 

3The mean of (J - K)p - (J - K)o is zero (-0.03 ± 0.29). 



The mean value of {Kp — K^) for the sys- 
tems with no IR excess is —0.28 ± 0.33. 
We find a similar offset for {Jp — Jo) of 
—0.30 ± 0.35^. We consider several pos- 
sibilities to explain the zeropoint offset. 
(1) Our predicted DA (V-K) colors are too red. 
Since we derive our DA K magnitudes from the 
predicted color {V — K)p and the published V 
magnitude, a difference in normalization in our 
synthetic V or K photometry of the DA models 
could cause the offset. However, the synthetic 
V and K mags and colors we derive for stars in 
the BPGS spectral library match their normal- 
izations and colors for each spectral type (see 
^ |4.1| above), and our {V — K) colors also match 



those in the spectral library of Pickles (1998) 
to within a mean difference of 0.01 ± O.lmag. 
(2) Our IR photometric zeropoint is too faint. 
Comparing the 2 DAs also observed by Zuck- 
erman & Becklin (1992), our K mags are O.lmag 
fainter in each case. While this difference is well 
within our errors, we cannot rule out a photomet- 
ric zeropoint shift of this order contributor to 
the mean zeropoint offset. (3) Unmodeled, large 
equivalent width IR absorption features in DAs. 
To account for a shift of 0.3mag in K, spectral ab- 
sorption features with equivalent widths of nearly 
1440A would be required that are not present 
in the DA models. Given that different absorp- 
tion features would be required in the J band to 
yield a similar zeropoint shift, we disregard this 
possibility. 

Even with this unexplained zeropoint shift, our 
determination of an IR excess, based on the sig- 
nificance of {Kp — Ko), are unaffected. However, 
our estimated companion magnitudes Kc and ab- 
solute magnitudes Mk^, may be up to 0.3mag too 
faint, leading in turn to spectral subclass esti- 
mates approximately half a subclass too late. 

The V magnitudes of the candidate compan- 
ions listed in Table ^ are in most cases signifi- 
cantly fainter than the visual magnitudes of the 
DA WDs. Since we elected not to observe in the 
IR those systems known to be composite by their 
optical spectra, this result is expected from our 
selection. 

4.4. Spatially Resolved Companions 

Finding charts for most of our sample have been 
published by Shara et al. (1993, 1997). To ac- 
count for visually resolved companions, we in- 
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Fig. 2. — Difference between predicted and observed magnitudes plotted against observed magnitude for IR-observed 
hot DA white dwarfs. Objects for which we determine a significant IR excess are shown as filled circles. 



spected images from the DSS (Shara et al. 1997), 
and also searched the USNOA-1.0 (Monet 1996) 
catalog to i? ~ 19.5 for all stars within 15" of the 
program star. Resolved stars within a maximum 
distance of 10" and within 3mag (in R) of the pri- 
mary are noted in Table as visual pairs. For the 
closest WD distances in our sample of ~ 20pc, 
the spatial resolution of ~ 2" in both the archival 
optical and new IR images, corresponds to 40 AU. 
On the other hand, a search limited to 10" corre- 
sponds to maximum separations of 200 AU. For 
the farthest WDs in our sample near 500pc, imag- 
ing is instead sensitive to resolved companions be- 
tween 1000 and 5000 AU. 

5. NOTES ON INDIVIDUAL OBJECTS 

REJ 0134-160: Based on radial velocity mea- 
surements (Schultz et al. 1996) and its low mass 
(NGS99, Finley et al. 1997) mass transfer is un- 
likely to have occurred in this binary DA-I-M star 
system. 

REJ 0148-253: Finley et al. (1997) claimed evi- 
dence for a companion based on H/5 emission and 
a slight red excess. Our result of an IR excess 
confirms the presence of the companion. 

REJ 0427+7406 This star appears to be a mul- 
tiple at near-IR wavelengths. Some nebulosity 
(diffuse emission) is also noted nearby. A find- 
ing chart is provided in Figure ^ 

REJ 0443-034: This is the hottest DA in our 
sample of IR-observed WDs. DA stars with such 
high temperatures should cool extremely rapidly 
and thus be detected only rarely. The suggestion 
by Finley et al. (1997) of a possible red excess is 
not borne out by our measurements. 

REJ 0512-004: As the least massive DA in the 
sample, this DA might be expected to have a com- 
panion, since very low mass WDs are thought to 
form in binaries (Marsh, Dhillon, & Duck 1995). 



The absence of IR excess indicates a very cool 
and/or degenerate companion at best. Our IR 
photometric method is insensitive to such objects, 
which are best detected through radial velocity 
variations (Saffer, Livio, & Yungelson 1998). We 
discuss the overall effect on our estimated binary 
fraction in section |^. 

REJ 0916-194: While its IR excess suggests a 
late-type companion, our estimate for its visual 
magnitude Vc is only ~ 2mag fainter than the to- 
tal V mag, so evidence for the companion should 
have been seen in the red end of its optical spec- 
trum. However, the only spectrum (from Vennes 
et al. 1997a) is of comparatively low S/N and 
extends only to 6040A. We expect that a higher 
S/N optical spectrum extending farther redward, 
or our own planned IR spectroscopy will readily 
reveal the companion. 

REJ 1036+460 (PG 1033+464): fsT photometry 
and ii"-excess from Zuckerman & Becklin (1992) 
are within O.lmag of the values we report here. 
However, their derived absolute red dwarf magni- 
tude Mk^ is fainter by 1.5mag. Rather than using 
detailed spectral fits, they had to use published 
color-estimated distance moduli that assumed 
uniform surface gravity for all white dwarfs, while 
in fact the absolute magnitude depends strongly 
on log (7. For the other object in common between 
our samples, REJ 1126+183 (PG 1123+189), we 
note an even more striking discrepancy (3.2mag in 
MkJ- They assume a distance of 30pc, while the 
best estimate from detailed fits of NGS is 135pc. 
Their resulting claim that the IMF is flat or in- 
creasing down to at least O.lMQshould be revis- 
ited with more accurate distance determinations. 

REJ 1043+490: This was suggested by 
Schwartz et al. (1995) possible multiple 

system based on a slightly red R — I color, but 
they suggested that further measurements were 
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needed to confirm this. The residuals of the DA 
model fit to the optical spectrum in NGS99 con- 
firmed a red excess for A > 4500A corresponding 
to a companion with V ~ 18.2 or My ~ 10.4 
at our distance of 370pc. However, our near-IR 
photometry is inconsistent with such a compan- 
ion, since a normal late M dwarf companion of 
this absolute magnitude and distance would yield 
J ~ 14.5 and K ~ 13.6. We note that the system 
is a visual binary, with a red companion 8" West 
of the DA. Contamination of the DA spectrum in 
NGS99 is thus the likely explanation. 

REJ 1235+233 (PG1232+238): The object is 
barely discernible in our K mosaic, so reported 
K photometry is poor. 

REJ 1711+664: This is a barely resolved visual 
pair, with a late-type star only ~ 2". 5 distant, 
which is included in our IR photometry (Mason 
et al. 1995). Jomaron (1997) also found the 
M-type secondary in residuals to the DA optical 
spectrum, classifying the former as type M3-M4, 
with no apparent H emission. Our IR photometry 
indicates a somewhat later type. 

REJ 2207+252: IR photometry for 2207+252 
was contaminated by a very bright, resolved ob- 
ject at 8.9". 

6. SELECTION BIASES 



Strong biases are likely to be inherent in any 
sample of EUV-detected WDs. The strongest 
is the selection against a high column of inter- 
vening interstellar material. Columns as low as 
5 X 10^^ cm~^ easily quench flux in the EUV be- 
low detection levels achieved in either the ROSAT 
WFC or the EUVE.'^ When compared to more 
distant areas in the Galactic plane, the local in- 
terstellar medium lacks cool dense material {n ~ 
0.1cm-3,T ~ lOOK). This 'local bubble' has a 
mean radius of 70 pc (Diamond et al. 1995). Ac- 
curate estimates of the intervening column to each 
star are difficult, requiring high S/N EUV spectral 
analysis well beyond the scope of this work. The 
global effect amounts to a direction-dependent 
ffux sensitivity limit (and therefore volume) that 
is strongly dependent on the fine structure of the 
ISM. 

As pointed out by NGS99, at a given temper- 
ature, WDs with lower masses are larger, conse- 
quently more luminous, and thus detectable to 
larger distances. For instance, at 30,000 K, a WD 
of O.5M0 is detectable at distances 30% larger 
than a O.7M0 WD. For many WDs at the max- 
imum distance sampled by most optical surveys, 
interstellar matter effectively absorbs all EUV ra- 
diation. Thus the dominant effect in EUV sur- 
veys is probably not a selection for massive WDs 



^For example, one optical depth is reached at column densities of 30, 5, 1, and 0.5 xlO^^ cm ^ respectively for the 
EUVE 100, 200, 400, and 6OOA detection bandpasses (Vennes et al. 1996). 
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al. 1997a, Finley, Koester & Basri 1997, NGS99). 
The 2MASS will provide (less sensitive, K < 14) 
IR photometry for brighter subsamples. 

WDs with masses less than the canonical core- 
helium ignition mass (O.49M0; Sweigart 1994) are 
likely to be in short-period binary systems, where 
their evolution was cut short by wind ejection dur- 
ing a common envelope phase (Marsh, Dhillon, & 
Duck 1995, Marsh 1995). Only one of the WD 
stars in our sample, REJ 0512-004, has a low mass 
(O.45M0), but it shows no signs of a companion. 
On the other hand, lower mass companions may 
be preferentially removed by merging, a process 
that is more effective the more massive the WD. 

It is thus of interest to check for any correlation 
between DA mass and IR excess. Figure § shows 
that no such correlation is detectable in our sam- 
ple, as confirmed by two correlation tests. Simi- 
larly, there is no correlation between the absolute 
K magnitudes of companions M^c and those of 
their corresponding DA white dwarfs.^ Two sam- 
ple statistical tests also reveal no significant dif- 
ference between the mass or absolute magnitude 
distributions of the DAs in the no-excess and ex- 
cess subsamples.^ Larger samples would put bet- 
ter constraints on these hypotheses. 

7. THE BINARY FRACTION 

Our WD sample demands that the initial pri- 
mary has already passed through the AGB phase, 
imposing a lower mass limit of about 2Mq on the 
WD progenitor. We thus probe portions of the 
(initial) mass ratio distribution that are quite dif- 
ficult to reach in unevolved systems. The binary 
frequency of high mass stars has proven difficult 
to determine for several reasons. Photometric de- 
tection of a faint low mass star next to a bright 
early type primary requires high precision pho- 
tometry. The use of radial velocity variations 
can also be difficult due to the large rotational 
speeds of many early-type stars - one reason why 
solar-type stars are preferred in current planetary 
searches. 

Out of our uniformly-selected sample of 57 
EUV-detected DA white dwarfs, including our 
new detections of IR excess, we find strong evi- 
dence for multiplicity in 20 systems (35%) within 

^Generalized Kendall's Tau and Cox Proportional Hazard model tests from the ASURV package (LaValley, Isobe & 
Feigelson 1992) show probabilities P of 18% and 72% respectively that a correlation is not present between ( Jp — Jo) and 
DA mass. A significant correlation requires P >99%. P of 26% and 36% respectively, are obtained for Mk^ vs. Mv^a- 

^Generalized Wilcoxon and logrank tests show that the two samples are inconsistent with being drawn from the same 
parent population at only the P ~ 75% level for DA mass and only the P ~ 93% level for My^A- 



as suggested by V97, but a selection against low- 
mass WDs due to a sample volume strongly af- 
fected by interstellar absorption. However, there 
is no direct bias induced by interstellar absorption 
on the observed binary fraction. 

There have been suggestions that detection of 
soft X-ray flux from hot white dwarfs may be 
biased against close binary systems. This could 
arise if heavy elements arising in a wind from 
the MS companion have accreted onto the surface 
of the white dwarf, enhancing atmospheric EUV 
and soft X-ray opacity (Barstow et al. 1993). 
By including only DA stars in our sample, we 
might thus be preferentially excluding some bi- 
nary systems. For instance, the hot white dwarf 
in the close binary RE 1016-053 is a DAO star 
whose mixed H/He composition is likely to be at- 
tributable to slow accretion from a red dwarf com- 
panion (Vennes et al. 1997b). However, Barstow 
et al. (1994) found that hot ROSAT (PSPC or 
EUV) detected WDs in non-interacting binaries 
with early-type stars showed no signs of abun- 
dance enhancements from accretion, their opacity 
sources following the same temperature pattern as 
seen for isolated WDs. Indeed, while those abun- 
dance enhancements may lead to significant X-ray 
opacity (e.g., Dupuis et al. 1995) in the hottest 
WDs, they remain strong EUV sources, particu- 
larly at A > 2OOA. 

Arguments for biases of opposite effect exist as 
well - that increased EUV flux is more likely from 
a WD that suffers reheating due to accretion of 
material from its companion. Since the DAs in 
our sample are not significantly variable, the in- 
teraction would have to be fairly weak or inter- 
mittent (e.g. via weak wind accretion). Further- 
more, Barstow et al. (1992) find for the DA-FK2V 
binary system V471 Tauri that the temperature 
change in such cases (with accretion rates of 
0.4 — 11 X 1O~^^M0 yr~^) is small and occurs only 
in the region near the hot spot of the accretion 
column. 

Measurement of these putative competing selec- 
tion effects would clearly benefit from larger sam- 
ples. A large number of EUV-detected DA white 
dwarfs with good constraints from atmospheric 
modeling now exist (Marsh et al. 1997, Vennes et 
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0,8 1 
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1.2 



Fig. 4. — J excess (difference between predicted and observed J magnitudes) plotted against DA white dwarf mass in 
solar units from the best-fit model atmospheres of NGS99. We see no evidence that IR excess correlates with DA mass. 



a 10" radius. Simple Poisson statistics on the 
subsample yields an error of 7%. We may also 
count multiple systems by assuming that any re- 
solved object within the 10" radius is at the 
same distance as the WD (i.e., a presumed com- 
panion), and thereby tally objects within a cho- 
sen limiting physical separation (projected onto 
the sky plane). Since the largest estimated 
distance to any of the WDs in our sample is 
470pc (for REJ 1235+233) we choose for com- 
pleteness a maximum projected physical sep- 
aration of lOOOAU. This removes only one sys- 
tem (1726+583), leaving the binary fraction es- 
sentially unchanged at 33 ± 7%. 

For comparison, among a complete sample of 
nearby {d < 5pc) M dwarfs, Leinert (1997) found 
26 ± 9% to have companions. In a meta-analysis, 
Fischer & Marcy (1992) found that 42 ± 9% of 
M dwarfs have main sequence companions. This 
is lower than the best incompleteness-corrected 
estimate for G dwarfs (~ 57%; Duquennoy and 
Mayor 1991), due in part to the smaller range of 
secondary companion masses available. 

As allowed by the faintness of WDs relative to 
main sequence companions, and the use of IR pho- 
tometry, the current study is sensitive to a wide 



range of companion masses, from late M to early 
A-type. The late-type companion masses corre- 
sponding to our estimates of absolute K magni- 
tudes range from about 0.5 down to O.IM©.^ 
Judging by known or estimated spectral types 
of unresolved companions and by derived abso- 
lute magnitudes of the presumed companions dis- 
cussed above, 90% of companions from Tables 
and 1^ are of spectral type K or later. 

We note that despite the significant improve- 
ment in sensitivity to low mass companions that 
our IR photometry affords, the binary frequency 
of our sample remains a lower limit for several 
reasons. First, very low mass companions can fall 
below our IR detection limit of ~ 16, with 
sensitivity decreasing with distance. Also, while 
we assume only a single degenerate in each sys- 
tem, double degenerates would not likely be dis- 
tinguished by our IR photometry. Saffer, Livio, 
& Yungelson (1998) derive a fraction of double 
degenerate systems that may be as high as 20% 
when they include an estimated correction factor 
of ~ 1.4 for selection efficiency in their sample. 

Some of the hot DAs with masses higher than 
expected from isolated evolution (~ I.IMq) have 
been hypothesized to represent a population of 



^Assumes a spatial resolution of 2" at a similar maximum distance of 500pc. Including only unresolved systems is 
intuitively appealing, but introduces a bias against companions in nearby systems. 

^Below about 0.2Mq, absolute magnitudes no longer correlate well with mass. There the relationship steepens consid- 
erably, in both V and K bands, and differs little as a function of metallicity (Baraffe et al. 1998; Henry et al. 1999). 
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coalesced double degenerate systems (e.g., Marsh 
et al. 1997). Such mergers would slightly reduce 
the observed binary fraction in our sample. While 
they might also help constrain the Type la super- 
nova formation rate (e.g.. Branch et al. 1995), 
the predicted peak of such a merger distribution 
should be near 0.9 Mq, which is not seen. 

Among our sample of hot DAs may be some 
close binaries that have passed through a phase 
of common envelope evolution. While these are 
sometimes called precataclysmic binaries, many 
will not become cataclysmic variables within a 
Hubble time (Ritter 1986). Based on morphologi- 
cal classification of planetary nebulae, and the as- 
sumption that axisymmetrical mass loss is caused 
by binary companions, Soker (1997) estimated 
that only 10% of PN had no interaction with 
any companion, stellar or substellar. Only 11% 
had interactions while avoiding a common enve- 
lope phase. For the case of close binaries, the 
observed distribution of companion masses may 
thus have been modified by evolution. Some sys- 
tems may have coalesced, while others may have 
grown through wind accretion, or by stable mass 
transfer. As the progenitors of novae, CVs, sym- 
biotics, Feige 24-type systems, and the intrigu- 
ing dwarf carbon (dC) stars (deKool & Green 
1995), these WD-I-MS systems warrant the im- 
proved constraints on their masses and abundance 
enhancements that IR spectroscopy of our sample 
will provide. For instance, the dC stars, which 
likely evolve into GH or Ba giants, are dwarfs 
whose strong carbon bands betray episodes of 
mass-transfer from the previously- AG B compan- 
ion, now the WD (Green 2000). There is recent 
evidence (Heber et al. 1993, Liebert et al. 1994) 
that dGs may show stronger abundance enhance- 
ments soon after the AGB phase, i.e., while the 
WD is still hot. IR spectra would test for this, 
and its potential ramifications for the depths and 
timcscalcs of mixing due to chemical potential 
gradients (e.g., Proffit & Michaud 1989). 

We are obtaining IR spectra of these systems at 
the KPNO 2.1m with GRSP to better constrain 
the mass and spectral type of the cool compan- 
ions presumably causing the IR excess. We await 
these results before pursuing more detailed dis- 
cussions of the mass function of the companions, 
the observed final and implied initial mass ratios 
of the systems, and detailed simulations of the 
effects of binary evolution. The binary fraction 
may be a function of the mass ratio, or of the pri- 



mary mass; these fundamental questions of stel- 
lar evolution remain open. If the observed binary 
fraction for DA WDs proves to be inconsistent 
with those of main sequence primaries in a simi- 
lar (initial) mass range, it could mean that age- 
dependent formation conditions (e.g., abundance, 
kinematics, or turbulence) strongly influence the 
binary fraction. 

Many thanks to Dick Joyce and the staff at 
KPNO for all their help during IRIM observ- 
ing runs. Thanks to Fred Ringwald for provid- 
ing a catalog of unresolved, hot-high-gravity with 
cool companion composites. The author grate- 
fully acknowledges support provided by NASA 
through Grant NAG5-1253, and Gontract NAS8- 
39073 (GXG). This research has made extensive 
use of SIMBAD. 
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APPENDIX 

PHOTOMETRY OF MOSAICED IMAGES 

For some targets, visual inspection failed to reveal the source on the individual images. We derive 
photometry these targets from a global mosaic of the field produced by shifting and combining individual 
frames. We use integer pixel shifts) and average the overlap regions between frames, accounting for 
background/sky variations by subtracting the median value of the frame from itself before averaging 
the fields. Integer pixel shifts do not necessarily conserve flux, but are employed here for two reasons: 
(1) the stellar images on our frames are sub-Nyquist sampled. And (2) the presence of significant 
dead (insensitive) space between pixels. The undersampled profiles of the star generally ensure that 
most of the stellar flux is contained within one pixel, and the presence of dead space combined with 
undersampled images makes it difficult to produce flux conserved stellar profiles in the first place. Thus 
we conclude that mosaicing images will not significantly alter our photometry. These assumptions 
are tested by comparing magnitudes derived from a mosaiced image to those from individual frames. 
Figures |A5| illustrated this comparison. Figure |A5| demonstrates that photometry from mosaiced images 
is consistent (within the errors) with photometry obtained from individual frames. 

The mosaiced images, however, only allow us to make an estimate of the formal (Poisson + detector 
noise) errors, and do not account for any systematic errors or any additional possible degradation of 
the S /N created by the photometry software. A comparison of the formal errors with those obtained by 
comparing the RMS magnitude variation from individual frame suggests that our errors are dominated 
by these systematic errors. The dashed line in Figure |A6| shows the magnitude of the formal errors. 



The solid dots are the observed RMS variations of magnitudes from individual frames. Clearly, formal 
errors underestimate the measured empirical estimates. The formal errors are derived using a modified 
form of the IRAF task ccdtime (modified to run on the IDL package) . To obtain a realistic estimate for 
the errors on mosaiced images, we modeled the effects of systematic errors by artificially increasing the 
readout noise of the detector and minimizing the relationship: 

i=l "obs 

where aobs are the observed errors, acai are errors calculated by increasing readout noise, and N is 
the number of stars in this experiment. The solid line in Figure shows an example of this fit. We 
use the best-fit "readout noise" to extrapolate the formal S/N to stars for which only mosaiced image 
photometry is available. 
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Fig. A5. — A comparison of the derived magnitudes in the field of PG 0904+511 using individual frames and a mosaic 
image, i^'-band magnitudes derived from an image mosaic and those from individual frames are compared for 5 stars in 
the PG 0904+511 field. The error bars show the 1-a deviation among magnitudes from individual frames. The dashed 
line shows the resulting relationship if the two quantities were identical. 
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Fig. A6. — The observed error (estimated as the 1-a deviation among magnitudes derived from individual frames) 

versus magnitude for the PG 0904+511 field. The dashed line shows the magnitude of the formal errors from S/N 
arguments. The solid line shows our error curve for the field (see text for details). These data suggest that we are 
dominated by systematic errors. 
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TABLE 1 

SAMPLE OF EUV-SELECTED DA WHITE DWARFS 



RE J Other Names V T^s log 5 M/Mq Refs Notes 





0007+331" 


GD2 


13 


85 


45653 


7 


77 





57 


1,2 






0134-160" 


GD984, PHL1043 


13 


96 


43722 


7 


70 





54 


1,2,9 


DA+dM 




0138+252 


PG0136+251 


15 


87 


38964 


9 


00 


1 


22 


1,1 


peculiar 




0148-253" 


GD 1401 


14 


69 


25707 


7 


70 





50 


1,2,6 


NS, DA+dM 





0235+034 


Feige 24 


12 


56 


62947 


7 


53 





54 


1,1,13 


DAe+dM 


0237-122" 


PHL 1400 


14 


92 


31570 


8 


38 





85 


1,2 




. 


0348-005 


GD 50 


14 


04 


38881 


8 


97 


1 


17 


1,2 




(N ■ 


0350+171 


V471 Tau 


13 


65 


34200 


8 


80 


1 


12 


1,1 


DA+KO 


^ \ 


0427+740" 




15 


9 


47549 


7 


84 





61 


3,2 


visual pair, 8" 


Ph. 

<\ 


0443-034" 




16 


9 


65140 


9 


12 


1 


29 


3,4 




0457-280" 


MCT 0455-2812 


13 


95 


50061 


7 


63 





53 


1,2 


NS 


: 


0512-004" 




14 


2 


31333 


7 


34 





45 


3,2 




0521-102" 




15 


81 


32775 


8 


53 





94 


1,2 




1 ■ 


0827+288 


PG 0824+289 


14 


22 


51934 


8 


00 





70 


1,1,12 


DA+dC+M3V 


> ■ 
: 


0841+032" 




14 


47 


37687 


7 


68 





52 


1,2 




0845+488 


HD74389B 


15 


5: 


39500 


8 


04 





65 


7,5 


DA+A2V 


^ ■ 


0902-040" 




13 


19 


22285 


7 


76 





51 


1,2 




^ ■ 


0907+505" 


PG 0904+511 


16 


54 


31791 


8 


04 





66 


1,2 




: 


0916-194" 




17 


3 


56400 


9 


12 


1 


29 


3,4 


visual pair, 6.4" 


.h/00 


0940+502" 


PG 0937+506 


16 





35511 


7 


62 





51 


3,2 




0957+852" 




15 


4 


50205 


8 


29 





82 


3,2 




1016-052 




14 


21 


53827 


8 


08 





65 


1,1,4 


DA+dM 


1 ■ 


1019-140" 




14 


93 


31100 


7 


85 





57 


1,2 




6 ■ 


1029+450" 


PG 1026+454 


16 


13 


35017 


7 


64 





51 


1,2 




in: 


1032+532" 




14 


45 


42785 


7 


89 





62 


1,2 




:as 


1033-114 


G 162-66, LTT3870 


13 


01 


23741 


7 


79 





52 


1,2 




1036+460" 


GD123 


14 


34 


28720 


7 


95 





61 


1,2,10 


DA+K 


> . 


1043+490" 




16 


23 


40461 


7 


87 





60 


1,2,14 


visual pair, 8" 


•i-H . 


1043+445" 


PG 1040+451 


16 


94 


48000 


7 


98 





63 


3,4 






1044+574" 


PG 1041+580 


14 


64 


29878 


7 


73 





52 


1,2 






1100+713" 


PG 1057+719 


14 


68 


40328 


7 


75 





55 


1,2 






1112+240" 


Ton 61 


15 


77 


39281 


7 


71 





53 


1,2 






1122+434" 


PG 1120+439 


15 


8 


26152 


8 


23 





75 


3,2 






1126+183" 


PG 1123+189 


14 


13 


52747 


7 


69 





56 


1,2,10 


DA+dM 




1128-025" 


PG 1125-026 


15 


73 


30227 


8 


15 





71 


1,2 






1148+183" 


PG 1145+188 


14 


33 


24717 


7 


84 





55 


1,2 






1235+233" 


PG 1232+238 


17 


4 


45639 


7 


77 





57 


3,2 






1236+475" 


PG 1234+482 


14 


38 


55570 


7 


57 





53 


1,1 


reclassified DB 




1257+220" 


GD153 


13 


38 


38324 


7 


71 





53 


2 






1316+295 


HZ 43 


12 


99 


49000 


7 


70 





56 


1,8 


DA+dM 




1336+694" 


PG 1335+701 


15 


4 


29067 


8 


27 





78 


3,2 






1426+500 




13 


6 


29800 


7 


98 





60 


3,4,6 


DA+MVe 




1431+370" 


GD 336 


15 


27 


33988 


7 


86 





58 


1,2 






1440+750" 




15 


32 


42400 


8 


54 





96 


3,4 






1446+632" 




16 


4 


37340 


7 


72 





53 


3,2 





TABLE 1— Continued 



REJ 


Other Names 


V 


^eff 


log 5 


M/Mq 


Refs 


Notes 


1629+780" 




13.03 


39882 


7.88 


0.61 


1,2,11 


DA+M4V 


1638+350'' 


PG 1636+351 


14.83 


34879 


7.90 


0.61 


1,2 




1643+411" 


PG 1642+414 


16.2 


28240 


8.15 


0.71 


3,2 




1650+403" 




15.83 


37698 


7.90 


0.61 


1,2 




1659+440" 


PG 1658+440 


14.62 


30510 


9.36 


1.31 


3,5 


magnetic DA 


1711+664" 




17.10 


47994 


8.78 


1.09 


1,2 


visual pair, 2.5" 


1726+583" 


PG 1725+586 


15.45 


52413 


8.14 


0.75 


1,2 


visual pair, 6.7" 


1800+683" 


KUV 18004+6836 


14.74 


43868 


7.74 


0.56 


1,2 




1820+580" 




13.95 


43296 


7.70 


0.54 


1,2 




1845+682" 


KUV 18453+6819 


15.5 


35599 


8.15 


0.72 


3,2 




2116+735" 


KUV 21168+7338 


15.0 


49777 


7.65 


0.54 


3,2 




2126+192 


IK Peg 


6.07 










DA+A8m 


2207+252 




14.58 


26129 


8.20 


0.74 


1,2 


visual pair, 8.9" 


2312+104" 


GD 246 


13.09 


51950 


7.76 


0.59 


1,2 





'^Ncw IR photometry presented in this work. 

NOTES: NS means not in official sample, but may have optical fits in NGS99 and/or IR data. 
Visual pairs within 10"and AR = 3mag are noted witli separations. 

REFERENCES: First reference is for V mag. 2d for Tea, logg, and M/Mq: 3d if hsted is for 
companion's spectral type and/or mass of previously known companion. If a multiple system has 
no third reference listed, it has several references available in SIMBAD. (1) - Marsh et al. 1997; 
(2) - NGS99; (3) - V from McCook and Sion; (4) - Vennes et al. 1997; (5) - Schmidt et al. 1992; 
(6) - FKB 1997; (7) - SIMBAD; (8) - Napiwotski, R. et al. 1993; (9) - Bues & Asian 1995; (10) - 
Green, Schmidt & Liebert 1986; (11) - Catalan et al. 1995; (12) - Green & Margon 1994; (13) - 
Vennes & Thorstensen 1994; (14) - Schwartz et al. 1995. 
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TABLE 2 
Summary of Observing Runs 





TTT Hntp 


kJlcllHJ.cliLl kJldlo 


kJldllQcllU. kJldlo 


V-'UllUl llUllo 






usea 


References^ 




KPOl 


04 Jan 96 


FS #6, #10, #19, LHS 


1,2 


Clear or light cirrus. 






254, LHS 325 






KP02 


05 Jan 96 


FS #10, #19, LHS 191, 


1,2 


Cirrus, but eventually 






LHS 2347 




very clear. 


KP03 


06 Jan 96 


FS #2, #10, #19, LHS 


1,2 


Heavy cirrus. 






2502 






KP04 


07 Jan 96 


FS #15, #10, #19, LHS 


1,2 


High thin cirrus. 






191 






KP05 


02 July 96 






Cloudy, no data. 


KP06 


03 July 96 






Cloudy, no data. 


KP07 


04 July 96 






Cloudy, no data. 


KP08 


05 July 96 






Cloudy, no data. 


KP09 


06 July 96 






Cloudy, no data. 


KPIO 


29 Oct 96 


FS #1,#6, #9, #12, #15, 


1 


Mostly clear. 






#29, #31 






KPll 


30 Oct 96 


FS #1, #10, #14, #29, 


1,2 


Photometric. 






GJ 1057, GJ 1002 






KP12 


20 Jan 98 






Cloudy. 


KP13 


21 Jan 98 






Fog. Telescope closed. 


KP14 


22 Jan 98 


FS #2, #6, #21, #24, 


1 


Clear. 






#18, #19 






KP15 


09 May 98 


FS #21, #23, #24, #27, 


1 


Clear. 






#29, #33 






KP16 


10 May 98 


FS #16, #20, #21, #23, 


1 


Mostly clear with moder- 






#35 




ate haze. 


KP17 


11 May 98 


FS #19, #21, #27 


1 


Hazy 



^FS=UKIRT faint star 

'^STANDARD STAR REFERENCES: (1) - CasaH & Hawarden 1992; (2) - Leggett 1992 
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TABLE 3 
Summary of Target Observations 





OViQPTviTi o" Tunnel 


1 'IL/ LCD 


0007+331 


KP02, KP04, KPIO 




0134-160 


KP02, KP04, 






KPIO 


jn3 


0148-253 


KP02, KP03, KP04 




0237-122 


KP04, 






KPIO 


kn9 


0348-005 


KPOl, KP02, KP03, KP04, 






KPll 


jn3, kn3 


0427+740 


KPOl, KP03, KP04, 






KP14 


kn9 


0443-034 


KP03, 






KPll 


jn9, kn9 


0457-280 


KP03, KPll 




0512-004 


KPOl, KP03, KP04, KPIO, 






KP14 


jn9, kn9 


0521-102 


KPOl, KP03, 






KPIO, 


jn9, kn9 




KP14 


kn9 


0827+284 


KPll 




0841+032 


KPOl, KPIO, 


jn9, kn9 




KP14 


kn9 


0902-040 


KPOl, KP04 






KPIO 


kn9 


0907+505 


KP03, KP04, KPll, KP14, 


jn9, kn9 




KP17 


knl8 


0916-194 


KP02 


jn9, kn9 


0940+502 


KPOl, 


jn9, kn9 




KP15, 


jcl, kcl 




KP16, 


jc2, jn9, kc2, kn9 




KP17 


kc9, knl8 


0957+852 


KP02, 






KP14, 


kn9 




KP16, 


jc2, jn9, kc2, kn9 




KP17 


kc9, kn9 


1019-140 


KP02 




1029+450 


KP02 




1032+532 


KPOl, KP04, 






KP14 


kn9 


1033-114 


KPOl, KP02, KP03, KP04 




1036+460 


KP02 




1040+451 


KP15, 


jcl, jnl2, kcl, knl2 




KP16, 


jc2, jn9, kc2, kn9 




KP17 


kcl2 
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TABLE 3 — Continued 



Target 


Observing Runs 


Notes" 


1043+490 


KP02, 






KP14, 


kn9 




KP15, 


jcl, jn9, kc2, kn9 


1044+574 


KP02 




1100+713 


KP02, 






KP17 


knl5 


1112+240 


KP16 


jc2, jn9, kc2, kn9 


1122+434 


KP02, 






KP16 


jc2, kc2, kn9 


1126+183 


KPOl, KP04 




1128-025 


KP02, 






KP16 


jcl, kcl 


1148+183 


KPOl, KP03 




1235+233 


KP14, 


kn9 




KP16, 


jc2, jn9, kc3, kn9 




KP17, 


kcl2, kn9 


1234+482 


KP14 


kn9 


1257+220 


KPOl, KP02, KP03 




1336+694 


KP03, KP04, 






KP14, 


kn9 




KP16 


jcl, jn9, kc2, kn9 


1431+370 


KP03, KP04 




1431+370 


KPIO, 


kn9 




KP16, 


kn9,jn9,jcl,kc2 




KP17 


kn9,kcl2 


1440+750 


KP14, 






KP16, 


jcl, jn9, kc2, kn9 




KP17 


kcl2, kn9 


1446+632 


KP14, 






KP16, 


jc2, jn9, kc3, knl5 




KP17 


kcl2, kn9 


1629+780 


KP15 


jcl, kcl 


1638+350 


KP15 


jcl, kcl, kn9 


1643+411 


KP15, 


jcl, jn9, kc2, kn9 




KP16, 


jc2, jn9, kc3, knl5 




KP17 


kcl2, kn9 


1650+403 


KP15, 


jcl, jn9, kcl, kn9 




KP16 


jc2, jn9, kc2, kn9 


1658+441 


KP14, 






KP15, 


jcl, jn9, kcl, kn9 




KP16, 


jc2, kn9, kc2, kn9 




KP17 


kcl2, kn9 


1711+664 


KP15 


jcl, jnl8, kcl, kn9 


1726+583 


KP15 


jcl, jn9, kcl, kn9 
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TABLE 3— Continued 



Target 


Observing Runs 


Notes" 


1800+683 


KP15, 


jcl, jn9, kcl, kn9 




KP16, 


jc2, jn9, kc2, kn9 




KP17, 


kcl2, kn9 


1820+580 


KPll, 






KP15 


jcl, jn9, kcl, kn9 


1845+682 


KPll, 






KP15 


jcl, jn9, kcl, kn9 


2116+735 


KPIO, KPll, 






KP15, 


jcl, jn9, kcl, kn9 




KP16, 


jc2, jn9, kc3, kn9 




KP17 


kcl2, kn9 


2127-221 


KPIO 


jn9, kn9 


2207+252 


KP15, 


jcl, jn9, kcl, kn9 




KP17 




2312+104 


KP02, KP03, KP04 




2353-243 


KPIO 


jn9, kn9 



'^The number of lOsec co— adds per exposure in both J and K bands is nominally 
6, and the number of dithered exposures is also 6, except as noted in the Notes. 

NOTES: 

XcY denotes that in the X bandpass the number of co-adds performed per ex- 
posure is Y (e.g., jcl means 1 co-add per exposure in J band). 

XnY denotes that in the X bandpass the number of combined dithered exposures 
is Y (e.g., kn9 means 9 dithered exposures combined in the iiTband). 



6 



TABLE 4 

Observed and Predicted Magnitudes and Colors of DA White Dwarf Stars 



RE Name 


J 


crj 


K UK {J-K) a^j_K) Mv {V-K)p Jp 


Kp 


Notes 








observed predicted 







0007+331 


14 


62 





02 


14.82 





07 


-0 


21 





07 


9 


06 


-0 


72 


14 


33 


14 


57 


0134 160 


12 


95 





01 


12.21 





01 





74 





01 


9 


00 


-0 


72 


14 


43 


14 


68 


0148-253 


12 


33 





01 


11.56 





04 





76 





04 


9 


98 


-0 


60 


15 


08 


15 


29 


0237-122 


15 


74 





06 


15.93 





12 


-0 


18 





14 


10 


68 


-0 


68 


15 


36 


15 


60 


0348-005 


14 


80 





06 


14.66 





52 





14 





52 


11 


46 


-0 


71 


14 


50 


14 


75 1 


0427+740 


15 


64 





03 


15.62 





12 





02 





12 


9 


13 


-0 


73 


16 


38 


16 


63 


0443-034 


17 


38 





05 


17.99 





07 


-0 


61 





08 


11 


23 


-0 


75 


16 


49 


16 


75 


0457-280 


14 


62 





02 


14.80 





03 


-0 


18 





03 


8 


69 


-0 


73 


14 


43 


14 


68 


0512-004 


14 


43 





03 


14.67 





01 





15 





40 


8 


81 


-0 


67 


14 


63 


14 


87 


0521-102 


16 


76 





20 


17.13 





13 


-0 


37 





24 


10 


87 


-0 


69 


16 


26 


16 


50 


0550-240 


16 


93 





10 


17.28 





41 


-0 


34 





42 


9 


03 


-0 


73 


17 


18 


17 


43 


0841+032 


15 


17 





02 


15.49 





01 


— 


33 





02 


9 


17 


-0 


70 


14 


93 


15 


18 


0902-040 


13 


70 





01 


13.90 





05 


-0 


20 





05 


10 


37 


-0 


53 


13 


54 


13 


72 


0907+505 


17 


16 





10 


17.54 





64 


-0 


39 





65 


10 


09 


-0 


68 


16 


98 


17 


22 


0916-194 


13 


94 





02 


13.64 





04 





36 





10 


11 


36 


-0 


74 


17 


78 


18 


04 


0940+502 


16 


73 





07 


16.97 





02 


-0 


24 





07 


9 


14 


-0 


70 


16 


45 


16 


70 


0957+852 


16 


49 





08 


16.72 





09 


-0 


23 





12 


9 


87 


-0 


73 


15 


88 


16 


13 


1019-140 


15 


17 





04 


15.45 





14 


-0 


28 





14 


9 


82 


-0 


67 


15 


36 


15 


60 


1029+450 


16 


92 





09 


17.54 





70 


-0 


62 





71 


9 


20 


-0 


69 


16 


58 


16 


82 


1032+532 


15 


38 





01 


15.45 





11 


-0 


07 





11 


9 


35 


-0 


72 


14 


92 


15 


17 


1033-114 


14 


80 





06 


14.66 





52 





14 





52 


10 


31 


-0 


56 


13 


38 


13 


57 


1036+460 


12 


62 





04 


11.89 





09 





72 





10 


9 


95 


-0 


65 


14 


76 


14 


99 


1043+490 


16 


86 





02 


16.95 





22 


-0 


10 





22 


8 


40 


-0 


72 


16 


71 


16 


95 


1043+445 


18 


00 





05 


17.48 





09 





52 





11 


9 


43 


-0 


73 


17 


42 


17 


67 


1044+574 


15 


32 





04 


15.73 





17 


-0 


42 





18 


9 


72 


-0 


66 


15 


06 


15 


30 


1100+713 


15 


51 





10 


15.78 





17 


-0 


27 





20 


9 


19 


-0 


71 


15 


14 


15 


39 


1112+240 


16 


37 





21 


16.51 





55 


-0 


14 





59 


9 


17 


-0 


71 


16 


24 


16 


48 


1122+434 


16 


12 





05 


16.39 





34 


-0 


27 





34 


10 


82 


-0 


61 


16 


20 


16 


41 


1126+183 


12 


74 





09 


12.04 





10 





70 





14 


8 


47 


-0 


73 


14 


61 


14 


86 


1128-025 


16 


18 





06 


16.48 





17 


-0 


30 





18 


10 


38 


-0 


67 


16 


16 


16 


40 


1148+183 


14 


96 





06 


15.22 





20 


-0 


26 





20 


10 


31 


-0 


58 


14 


71 


14 


91 
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TABLE 4— Continued 



KrL iName 


■ 

J 


r 




K cjK {J- K) 
observed 




-K) 


My 


{V - K)p Jp 
predicted 


TV' 
K 


p 


Notes 


1235+233 


18 


31 





11 








-2 


36 





28 


9 


06 


-0.72 


17 


88 


18 


12 


2 


1236+475 


14 


91 





06 


14.90 





16 





01 





17 


8 


37 


-0.73 


14 


86 


15 


11 




1257+220 


14 


17 





06 


14.43 





02 


-0 


26 





05 


9 


21 


-0.71 


13 


84 


14 


09 




1336+694 


15 


85 





02 


16.12 





04 


-0 


28 





04 


10 


66 


-0.65 


15 


82 


16 


05 




1431+370 


15 


97 





04 


16.37 





16 


-0 


40 





17 


9 


65 


-0.69 


15 


72 


15 


97 




1440+750 


17 


21 





11 


16.64 





12 





57 





16 


11 


39 


-0.70 


15 


77 


16 


02 




1446+632 


16 


77 





11 


17.77 





07 


-1 


00 





12 


9 


26 


-0.70 


16 


86 


17 


10 




1629+780 


10 


96 





03 


11.19 





04 


-0 


23 





05 


9 


42 


-0.71 


13 


49 


13 


74 




1638+350 


15 


55 





11 


15.57 





64 


-0 


01 





64 


9 


67 


-0.70 


15 


28 


15 


53 




1643+411 


16 


86 





10 


17.46 





05 


-0 


61 





11 


10 


53 


-0.64 


16 


62 


16 


84 




1650+403 


16 


60 





08 


16.86 


1 


05 


-0 


26 


1 


05 


9 


54 


-0.71 


16 


29 


16 


54 




1658+441 


15 


47 





10 


15.68 





02 


-0 


21 





10 


12 


70 


-0.67 


15 


06 


15 


29 




1711+664 


15 


02 





01 


14.19 





18 





83 





18 


10 


84 


-0.73 


17 


58 


17 


83 




1726+583 


16 


40 





26 


16.17 


1 


51 





23 


1 


53 


9 


56 


-0.73 


15 


93 


16 


18 




1800+683 


15 


54 





04 


15.87 





04 


-0 


33 





06 


9 


06 


-0.72 


15 


21 


15 


46 




1820+580 


14 


07 





04 


13.42 





01 





65 





04 


9 


01 


-0.72 


14 


42 


14 


67 




1845+682 


14 


85 





10 


14.37 





19 





47 





22 


10 


06 


-0.70 


15 


95 


16 


20 




2116+735 


15 


54 





01 


15.93 





06 


-0 


39 





06 


8 


74 


-0.73 


15 


48 


15 


73 




2207+252 


14 


34 





08 


14.06 





22 





28 





23 


10 


77 


-0.61 


14 


98 


15 


19 


3 


2312+104 


13 


81 





04 


14.06 





01 


-0 


14 





75 


8 


88 


-0.73 


13 


57 


13 


82 





NOTES: Objects with significant IR excess, are shown in boldface. 
-^Faint IR standard star - photometry from CasaH & Hawarden (1992) 
^Object is barely discernible on K band image; K photometry unreliable. 
^IR photometry known to include nearby optically resolved object at 2.2" 



8 



TABLE 5 

Derived Properties of Candidate Cool Companions 



REJ Name 


Jc 






M Subtype'" 




Notes 










(YD) 


(OD) 


(OD) 




0134-160 


13.3 


12.3 


7.4 


4.5 


3.5 


17.4 


1 


0148-253 


12.4 


11.6 


6.9 


4 


3 


16.4 


1 


0427+740 


16.4 


16.2 


9.4 


6.5 


6 


23.5 


3 


0916-194 


14.0 


13.7 


7.7 


5 


4 


19.2 




1036+460 


12.8 


12.0 


7.6 


5 


4 


17.3 


1 


1126+183 


13.0 


12.1 


6.5 


3.5 


2.5 


16.6 


1 


1629+780 


11.1 


11.3 


7.7 


5 


4 


16.6 


1 


1711+664 


15.1 


14.2 


8.0 


5.5 


5 


20.0 


2 


1820+580 


15.5 


13.8 


8.9 


6 


5.5 


21.0 




1845+682 


15.3 


14.6 


9.2 


6.5 


6 


21.9 





Approximate spectral type from Mk^ and Leggett (1992) Table 6, for 
Young Disk and Old Disk stars. 

''Approximate V magnitude of companion assuming Old Disk popula- 
tion, from Mk, and {V - K) of Leggett (1992) Table 6. 

NOTES: 1) known binary; see Table 1. 

2) IR photometry known to include nearby marginally optically-resolved 
object at separation indicated in Table 1. 

3) This star is spatially resolved at near-IR wavelengths, with a compan- 
ion within < 2". Some nebulosity (diffuse emission) is also noted nearby. 
The more distant companion noted in Table 1 is not included in the IR 
photometry. 
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